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We shortly review recent successes in applying Bose- Einstein interferometry in heavy ion colli- 
sions and the proceed to some model calculations for 3-dimensional Bose- Einstein correlation 
functions in e + e~ collisions at the Z° pole. 



1 Theoretical Overview 

Bose-Einstein correlations (BEC) are a phase-space phenomenon: Symmetrization of the mul- 
tiparticle wave function affects the measured n-particle coincidence spectra and leads to an 
enhancement relative to the corresponding product of independent 1-particle spectra, if the 
emitted particles are close in phase-space (i.e. they occupy the same elementary phase-space 
cell). The spatial length of the elementary phase-space cells is limited by the geometric size of 
the source of particles with the considered momentum. The larger this size, the narrower these 
cells are in momentum space. By tuning the relative momenta and watching the onset of BEC 
effects one can thus measure the spatial length of the elementary phase-space cells and thereby 
the size of the source. 

Wigner Functions. A description of BEC effects among n particles thus involves the n- 
particle phase-space density. Since we are discussing a quantum mechanical phenomenon, we 
are not talking about a classical phase-space density (which has directly a probabilistic interpre- 
tation), but about the Wigner density (which is positive definite only when averaged over many 
elementary phase-space cells). If the particles are emitted independently, the (unsymmetrized) 
n-particle Wigner density factorizes, and all n-particle coincidence cross sections are expressible 
through the single-particle Wigner function S(x,p). The assumption of independent particle 
emission is justifiable in heavy ion collisions where the many unobserved particles serve as a 
reservoir for all kinds of conserved quantities. In e + e~ collisions this is much less obvious and 
needs to be tested experimentally. 



Correlation Function. As long as the source has sufficiently low phase-space density that 
multi-particle symmetrization effects are dominated by two-particle exchange terms, the two- 
particle correlation function C(q,K), defined as the ratio of the 2-particle coincidence spec- 
trum P2(p a ,Pb) and the product of single-particle spectra P\{p a )Pi{pb) with q = p a —Pb and 
K = (Pa+Pfc)/2, is given bylT 
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Here J x = J d 4 x, q° = E a -E b , K° = (E a +E b )/2, and 
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The normalization Af depends on the multiplicity distribution via@ J\f = (n(n—l))/(n) 2 . In 
heavy ion collisions usually J\f ~ 1 . Due to the mass-shell constraint^ q° = (3 • q (where 
(3 = K /K° faK/ Ek is the velocity of the particle pair) the Fourier transform in ([!]) is not 
invertible: the separation of temporal and spatial aspects of the emission function S(x, K) re- 
quires additional model assumptions which must be provided by a physical picture of the time 
evolution of the source until freeze-outE 

The Reduced Correlator. While (jlj) goes to 2J\f at q = 0, real correlation functions usually 
approach a smaller value jV(l+A) with X(K) < 1. Possible reasons are partial phase coherence in 
the source and decay contributions from long-lived resonancesB To account for this one rewrites 
© as 
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The reduced correlator lC TC &{q,K) is given by the last term in (|l]) which contains the informa- 
tion about the space-time structure of S(x,K). To isolate it one constructs C(q,K) from the 
measured 1- and 2-particle cross sections, applying the Coulomb correction, determines A/ - and 
X(K) from the limits q — > and q — > oo, divides by M and subtracts the 1, and finally divides the 
result by X(K) and the measured ratio of single particle cross sections Pi(K) 2 /Pi(p a )Pi(pf ) ). 
For large sources like those in heavy ion collisions this ratio is close to unityp but for small 
sources like those in e + e~ it can contribute significantly to the g-dependence of C(q,K); it is 
then important to divide it out before trying to extract the source size. So far we have seen no 
data analysis where this is done! Instead, one usually extracts the size directly from K.(q,K), 
without dividing out the 1-particle spectra. As we will see, this can be quite misleading. 

Source Radii from BEC. One usually characterizes^ the source function S(x, K) by its norm, 
center and space-time variances (widths) , all of which are generally functions of the momentum 
K of the emitted particles. In this "Gaussian approximation" the reduced correlator reads 
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are the space-time variances of the emission function (effective source sizes). Different conven- 
tions for resolving the mass-shell constraint q° = (3 ■ q and expressing (||) in terms of three 

"We here neglect Coulomb final state interactions since methods are known to correct the data for themlll 



indepenent components of q lead to different Gaussian parametrizations for the correlatorcl The 
corresponding Gaussian width parameters, the "HBT (Hanbury Brown - Twiss) radii", are then 
combinations of the variances (XuX u )(K) and thus functions of the pair momentum K. 

2 Bose-Einstein Correlations in Heavy Ion Collisions 



Pb+Pb collisions at the SPS it was found that the pion emitting source is a rapidly expanding 
fireball in approximate local thermal equilibrium which at decoupling has a temperature of 
about 100 MeV and expands nearly boost-invariantly in the longitudinal direction while the 
average transverse expansion velocity is a bit larger than half the light velocity. The collective 
expansion manifests itself in a strong and characteristic dependence of the space-time variances 
{x^x v ) of the effective source S(x, K) on the pair momentum K. This implies a corresponding 
.ff -dependence of the HBT radii extracted from (Q). The pion emission process lasts only for 
about 2-3 fm/c but it doesn't begin until at least 6-8 fm/c after the collision. Freeze-out thus 
is a rather sudden process at the end of an extended rescattering and expansion stage. It is 
important to stress that the separation of longitudinal and transverse flow and access to the 
emission duration (P) is only possible in a full-fledged 3-dimensional and /^-dependent analysis 
of the correlation function C(q,K). Projections to lower dimensionality (e.g on qf nv ) lead to 
uncontrollable and unrecoverable loss of information. 

3 Bose-Einstein Correlations in e + e Collisions 

As stated in Sec. 1, to compute Bose-Einstein correlations one needs information on the Wigner 
phase-space density of the source. Going simulation programs of particle production in high- 
energy e + e~ collisions like PYTHIA, JETSET and HERWIG provide only momentum-space 
information on the produced particles. This is not enough to calculate BEC effects. Different 
methods have been suggested to provide the missing coordinate-space information, either di- 
rectly or indirectlyB We previously studiedEl BEC in VNI which studies the time evolution of 
the collision in phase-space. Here we present some very early results based on a phase-space 
versionlil of JETSET 7.4 which provides both the momenta and production coordinates for the 
produced particles. Our version of this code distributes the transverse distance of the produc- 
tion points from the central string axis according to a Gaussian with rms radius of 0.78 fm while 
S. Todorovova's version0 puts the production points right on the string axis. This latter proce- 
dure is inconsistent with the uncertainty relation, and we found accordingly^ that it produces 
correlation functions which rise as a function of q s instead of decaying. 

The algorithm for computing the correlation function from the positions and momenta of the 
generated pions is described elsewherepl we use the "classical" algorithm without wave packet 
smearingB In order to test the space-time structure of the events generated by JETSET and 
the BEC afterburner, we begin with a simple event topology (e + e~ — > Z° — > qq — > 2 jets) and 
consider only directly produced pions, thus avoiding the multiscale problems associated with 
longlived resonance decays. We analyse the correlation function in a Cartesian coordinate system 
where the longitudinal (I or L) axis is along the direction defined by the relative momentum of 
the initial qq pair (~ jet axis), the outward (o or T) direction is defined by the transverse pair 
momentum Kt, and the sideward (s) axis points in the third direction. 

The top two left panels of Fig.|l] show the correlator in the side direction. The reduced 
correlator /C re( j is seen to be independent of and always reproduces the input rms width 
of the string: R s = r vms /\/2 = 0.55 fm. In contrast, /C does depend on K^, and for small Kt 
it produces smaller HBT radii (0.31, 0.41, 0.46 and 0.50 fm at K T = 0, 0.3, 0.5 and 1.0 GeV, 
respectively). This effect is an artifact induced by the ratio of 1-particle spectra in (||); it 




Due to space reasons we will be very 



short - detailed discussions can be found elsewhere crEi For 





matters since the real radius is so small, producing significant errors if not divided out. For R Q 
and Ri, on the other hand, its effect is in our calculation nearly negligible: these radii come out 
much larger than R s . This, however, points to another problem: longitudinal HBT radii of up 
to 5 fm are incompatible with the data which give only about 1 fm (see the experimental talks 
in this session)! The problem seems to be connected with the large emission time duration At 
of up to 3 fm/c at low K^. This parameter, which reflects the proper time distribution of string 
breaking processes in JETSET, is not fixed by 1-particle spectra, but it is seen to seriously 
affect the 2-particle correlations. We are presently trying to fix this problem. At this moment 
we can only say that the version of JETSET used by us disagrees with experiment at the level 
of 2-particle correlations. 
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Figure 1: First two panels: The correlation function in side-direction (q = qi=0) for pairs with Kl = (such that 
the Roi cross term vanishesa) and Kt=0 and 0.5 GeV, respectively. Third panel: The emission time duration 
as a function of Kt for Kl, = 0. Second row: R a , R , and Ri as functions of Kt for Kl =0. 
We checked that the HBT radii correctly reproduce the rms widths of the space-time scatter plots of the produced 

pions in the appropriate K- windows. 
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